Nature Physics

Cell density is a fundamental biophysical property of living cells, yet the dynamical architecture underlying density homeostasis remains unknown. Here, we reconstruct the geometry, dynamics, and control structure governing cellular density regulation from single-cell mass-volume measurements acquired during recovery from osmotic perturbation. We identify a persistent anisotropic attractor in mass-volume state space that permits coordinated density-preserving fluctuations while strongly suppressing density-changing deviations. Optimal transport analysis reveals directed recovery trajectories confined to this attractor, enabling inference of an effective density-restoring control law. Density deviations are associated with coordinated changes in biomass accumulation and volume expansion that together generate an effective density-restoring control law. Moreover, density error contains significant predictive information regarding future corrective behavior, revealing an information-rich component of cellular regulation. Together, these findings establish cellular density homeostasis as an emergent dynamical control process and demonstrate how attractor reconstruction, optimal transport, and information theory can be combined to uncover the organizing principles of biological homeostasis.

Neuronal dendrites form highly branched morphologies supported by spatially patterned microtubule arrays that are essential for intracellular transport and structural stability, yet how this organization is established and maintained remains unclear. Here we analyze microtubule behavior in Drosophila sensory neurons using live-cell imaging, deep learning-based analysis, and genetic perturbations. We find that microtubules are shorter and fewer in thin distal branches and progressively longer and more abundant in thicker proximal branches. Locally, polymerization slows in narrow branches due to soluble tubulin depletion, while catastrophe frequency increases near branch points upon cortical encounter. A physical model incorporating these morphology-dependent constraints reproduces both local microtubule dynamics and global organization. Inverting the textbook view that microtubule dynamics shape cell morphology, our results demonstrate that, on short timescales, cell morphology instead constrains microtubule dynamics, thereby determining large-scale cytoskeletal organization.

How mesoscale cargos such as biomolecular condensates are transported and spatially positioned within cells remains a fundamental question in cell biology. During Caenorhabditis elegans zygote polarization, a substantial fraction of P granules undergoes persistent anterior-to-posterior motion that cannot be fully explained by dissolution-recondensation. Here, combining microfluidic reconstitution, quantitative micropipette aspiration and live-cell imaging, we show that a MEG-3 concentration gradient generates interfacial-tension asymmetry and drives the directed transport of PGL-3-RNA co-condensates in vitro and P granules in vivo through the Marangoni effect. This framework enables parameter-free quantitative predictions of condensate transport. Orthogonal molecular perturbations abolish or reverse directional migration, directly linking interfacial properties to transport behavior and establishing the specificity of this mechanism. Together, our findings identify the Marangoni effect as a previously unrecognized mechanism for intracellular transport of biomolecular condensates and provide direct causal evidence that condensate emergent properties can contribute to physiological organization.

The emergence of treatment resistant mutants is one of the primary contributors to therapy failure in cancer. Evolution-based adaptive therapies aim to exploit the competition between sensitive and resistant cells to prolong therapeutic control. However, in the densely packed tissues of solid tumors, competition is not only determined by cell-intrinsic factors but also a result of local physical interactions between neighboring cells. The successful development of physics-informed adaptive therapies therefore hinges upon the bottom-up integration of tissue mechanics and experimentally grounded models of tumor evolution. Here, we show how population-level treatment response emerges from mechanical interactions between individual cells. Combining genetically tailored therapy-mimicry experiments in spatially structured 2D melanoma colonies with mechanistic, particle-based tissue simulations as digital twins in a real-to-sim-to-real framework, we reveal how treatment-induced changes in mechanical crowding shape a cascade of resistance expansion, spatial competitive escape, and consequent therapy failure. We then link these underlying physical mechanisms to population-level observables in an abstracted mathematical model for reinforcement-learning-based therapy optimization and a distilled treatment-decision-heuristic, which we test experimentally. Together, our results identify mechanical crowding and escape as governing principles of resistance-mediated therapy failure in solid tumours and provide a basis for mechanistically informed, evolution-based treatment strategies.

The actin cytoskeleton drives shape changes and transport in cells and supports mechanical signal transmission. However, how cells control and use active cytoskeletal flows at the mesoscale is not well understood. We reconstituted an active cytoskeleton in water-in-oil emulsion droplets of Xenopus laevis egg extract and observed, above a critical droplet size, the emergence of a 3D radially convergent steady-state flow of polymeric actin, maintained by continuous actin turnover. The flow condensed lipid-rich cellular debris into a centered inclusion. Steady-state F-actin density and flow velocity profiles roughly collapse onto scale-invariant master curves. This behavior can be explained by a physical model representing the network as an isotropic active viscous fluid with a percolation threshold. The contracting network behaves as an active swimmer with complex internal dynamics that centers itself and the central inclusion inside the droplets without physical boundary attachment. Active contraction, crosslinking and polymerization dynamics in an actin network can thus generate cell-scale flow patterns that sense the confining geometry and external signals and exert forces that are likely sufficient to move and localize organelles in cells.

Vascular remodeling responds to local wall shear stress, but transport efficiency belongs to the network as a whole. We connect these scales by minimizing dissipation under a fixed maintenance budget, with vessel cost proportional to mass raised to an exponent b. At an optimum, dissipated power is proportional to vessel cost. This yields the shear set point |{tau} | {propto} r^(b-1){ell}^(b-1)/2 for radius r and length {ell}, with uniform shear recovered at b = 1. The associated local adaptation rule decreases the global functional, although the topology it reaches retains a memory of development. Simulations with 400 and 700 sinks show that staged growth improves transport mainly by allowing pruned connections to regrow. Fluctuating demand sustains loops without increasing the maintenance budget. The local condition fixes allocation on a given support, whereas growth and variable demand shape the topology that survives.

Collective migration is often attributed to complex cell-cell communication, yet coordinated multicellular dynamics can also emerge in organisms lacking social interactions. Continuum descriptions successfully capture such behaviors but generally rely on phenomenological transport coefficients whose microscopic origin remains elusive. Here we show that these coefficients emerge from cell-cell collisions and can be independently determined from experiments. We investigate the aerotactic dispersal of confined colonies of the solitary amoeba Acanthamoeba castellanii, in which oxygen consumption continuously reshapes the local environment to generate self-organized guidance cues. Combining quantitative videomicroscopy, oxygen mapping with phosphorescent sensing films, and oxygen-controlled microfluidics, we directly measure both the oxygen landscape and the aerotactic response governing colony expansion. We then develop a Patlak-Keller-Segel (PKS) -type drift-diffusion model whose transport coefficients are entirely determined from independent experiments and microscopic-to-macroscopic modeling approach, without adjustable parameters. The model quantitatively reproduces the observed dynamics, including the coexistence of a fast ballistic front and a trailing transient accumulation band. We show that both diffusion and aerotactic transport are modulated by the same collision-controlled mechanism through density-dependent changes in cell persistence. More generally, our results establish transport coefficients as measurable emergent properties bridging microscopic cell dynamics and macroscopic collective behavior, providing a predictive framework for oxygen-guided migration in living active matter.

A foundational tenet of modern biophysics dictates that protein-solvent interactions and charge-transport dynamics operate strictly within microscopic, sub-nanometer boundaries where long-range coordination is natively suppressed by rapid room-temperature thermal dephasing. Here, we report an unexpected phenomenon that fundamentally challenges this paradigm: the long-distance, quantum-coherent communication between isolated protein samples in aqueous solution across macroscopic spatial separations. By mapping the distance-dependence of a weak, protein-emitted electromagnetic field on the response dynamics of an adjacent complex, we demonstrate that this inter-protein communication establishes a stable, intermolecular vibronic coupling operating across millimeter-scale distances up to at least 40 mm. This field interference induces a robust quantum coherence that lowers the local potential energy surface, transforming the coupled protein-solvent network into an un-attenuated, low-barrier energy transfer conduit. This stabilized pathway launches a topologically protected, non-dispersive soliton that exhibits exceptional phase stability against room-temperature thermal dephasing. These findings provide definitive empirical evidence of a sophisticated, field-driven quantum infrastructure embedded within a biological macromolecule. More explicitly, they reveal that a protein can possess an inherent capacity to preserve coherent wave propagation over macroscopic scales previously deemed impossible in room-temperature aqueous environments.

The actin cytoskeleton is a self-assembled material that acts as the mechanical scaffold of the cell. Within the cell, protein cross-linkers organize actin filaments into assemblies with various microstructures that support distinct physiological functions. These protein cross-linkers differ in structural features. To decouple the role of individual cross-linker properties in actin assembly microstructure and mechanics, we synthesize programmable cross-linkers using a DNA linker coupled to binding domains. We find that high affinity cross-linkers induce actin bundling, while low affinity cross-linkers form actin networks. Myosin II motors remodel and contract all cross-linked networks tested in distinct ways from uncross-linked networks. We find that networks with high affinity cross-linkers have different myosin induced dynamics than networks with low affinity cross-linkers. This research expands our understanding of how cells regulate their mechanics and microstructures, while also providing insight towards the structure-property relationship in complex soft materials.

The back-and-forth (bidirectional) motion of Endosomes and Phagosomes along microtubules is fundamentally important for nutrient absorption, receptor recycling, pathogen degradation and many other downstream cellular pathways. Bidirectional motion requires opposing force from the kinesin and dynein motors, causing intermittent Tug-of-Wars (ToWs) followed by runs where winning motors drive their cargo until the next ToW. Dynein displays a counter-intuitive behaviour of catch bonding tightly to the microtubule against opposing force, but why such a response is needed remains unknown. Here we force dynein to catch-bond using an optical trap when dynein-teams are transporting phagosomes in-vitro and inside cells. An unexpected form of collective inactivity emerges, where larger dynein-teams remain inactive for longer after they have been catch-bonded. Kinesin exhibits a high probability of reversing the cargo during such periods of inactivity of dyneins. Therefore, the inactivity of dyneins allows cells to auto-regulate the effect of opposing molecular forces for the purpose of achieving bidirectional motion of endosomes/phagosomes. We propose that an inactivating conformational change is induced when dynein detaches against strain from its catch-bonded state. Collective inactivity follows because each inactive dynein exerts an entropic force on the other to delay the entire team from being activated. We capture the dynamics of this process through a statistical model of Entropically Constrained Reactivation, obtaining predictions in excellent agreement with experimental data. The emergent behaviour observed here possibly prevents a runaway dynein situation inside cells, its simple entropic origin suggesting broader relevance to catch-bonded receptor ligand interactions operating under conditions of cellular crowding.

Multi-cellular pattern formation is a key process in organism development, regeneration, and disease. Such multi-cellular patterns include, for instance, cell density and gene expression patterns. Here, we develop a cell-track-based method that quantifies how much individual cellular processes contribute to the emergence of large-scale patterns. These cellular processes include changes in cellular expression levels, cell motion, cell divisions, and apoptosis. We validate our method using synthetic data, before demonstrating it on two sets of experimental data. Our method is purely geometric and does not require any mechanistic assumptions about the underlying pattern-forming processes, making it very generic and robust. It thus enables unbiased, data-driven study of the mechanisms driving pattern formation across a broad class of systems.

Lumens populate most of our organs and their size, which dynamically changes throughout the day and lifetime, can influence their function. During mammalian development, the first lumen appears in the blastocyst. Surface epithelial cells pump ions to build an osmotic gradient that draws water, which inflates the embryo and stretches the epithelium as well as the zona pellucida, an elastic shell protecting the embryo. Blastocyst expansion is not steady and instead undergoes cycles of expansion interrupted by collapses during which the lumen leaks out. Early studies identified the critical role of osmotic pressure difference in regulating the expansion of the blastocyst whereas recent studies proposed that mechanical stresses would set the size of the embryo. However, the relative contributions of osmotic and mechanical stresses to lumen expansion dynamics remain unclear. Here, we investigate this using microfluidics, biophysics, machine learning and modelling on mouse embryos. Breaking down the expansion dynamics of mouse embryos, we discovered that fluid pumping rate accelerates, compensating for the increased frequency of collapses interrupting episodes of expansion. We reasoned that fluid transport could be rate limited by osmotic pumping, water permeability and/or mechanical constraints. Removing the mechanical constraint from the zona pellucida only marginally affected these expansion dynamics. Forcing fluid exchange with osmotic shocks we measured that water permeability is large enough not to restrict fluid accumulation rates. Finally, we estimate ion pumping for mouse blastocysts and interfere with it using inhibitors to measure how it quantitatively sets blastocyst expansion rates. Together, we propose that osmotic pressure differences, regulated by fluid transport properties, dominate over mechanical stresses in controlling blastocyst expansion dynamics.

Myelin membranes are tightly connected to the axon, forming an electrical route for rapid saltatory conduction of action potentials from node to node. Beneath the myelin, adjacent to the node, lie voltage-gated potassium (K+) channels, whose fundamental role in conduction remains poorly understood. Constrained by biological data, we developed a finite element model of a node of Ranvier encompassing an adaxonal membrane separating oligodendrocyte cytoplasm from the periaxonal space, and solved the full coupled set of physical equations governing electrodiffusion and voltage profiles at the nanoscale. Remarkably, we find that nodal action potentials invert at the inner adaxonal myelin membrane. The hyperpolarization-driven opening of myelinic Kir4.1 channels produces an instantaneous, powerful K+ prebuffering, clamping the periaxonal K+ concentration and enabling high-frequency action potential generation. Together, these findings indicate that the periaxon maintains an electrical feedback loop with the node, and suggest a key evolutionary advantage of juxtaparanodal K+ channel positioning.

Extrachromosomal DNAs (ecDNAs) form highly heterogeneous condensates in cancer cells that drive oncogene overexpression, yet how structural variability coexists with stable gene regulation remains unclear. Here, we develop a minimal polymer physics model of MYC-harbouring COLO320-DM ecDNAs, where BRD4-like complexes bind and bridge cognate sites along ecDNA rings. Above a critical binder concentration, ecDNAs phase separate into condensates exhibiting diverse conformations because of their thermodynamic folding degeneracy. Despite this variability, condensates retain conserved interaction scaffolds that give rise to reproducible contact patterns, including in-trans associated domains (I-TADs), genomic regions enriched in intermolecular regulatory contacts between distinct ecDNAs. We find that condensate 3D architecture follows universal scaling relations with ecDNA copy number, n, remaining robust to model parameter changes. Regulatory contacts within I TADs increase linearly with n, yet they are one order of magnitude stronger than in size matched control regions outside I TADs, whereas their relative fluctuations are markedly suppressed as n increases. This scaling produces enhanced, low-noise regulatory environments for oncogenes embedded within I-TADs, such as PVT1-MYC fusions, whereas the canonical MYC copy, located outside, is less amplified as experimentally observed. Our findings reveal universal polymer physics principles underlying ecDNA condensate organization, offering a mechanistic basis for selective oncogene amplification and potential advantages in cancer progression.

Influenza virus attaches to cells by binding to sialylated glycans. Subsequent endocytosis of the virus and endosomal maturation leads to a conformational change in the hemagglutinin coat protein that ultimately results in membrane fusion and release of the viral genome into the host cytoplasm, where it is trafficked to the nucleus and commences replication. Individual hemagglutinin-glycan binding affinities are weak (Kd approximately mM), so multivalency is key to viral attachment. Receptor density or clustering has previously been shown to affect influenza binding, but the overall spatial heterogeneity of both plasma and endosomal membranes has not been fully addressed. Here, we show that nanoscale lipid clustering can increase local concentration of glycosphingolipids and greatly enhance influenza viral binding in a way that macroscopic phase separation does not, and that binding avidity is dependent on the spatiotemporal scales of this clustering. This clustering is driven by sphingolipids and cholesterol and the nanoscale lipid environment of these clusters persists to affect the viral entry process. This finding demonstrates how fine receptor patterning can drive infection and illustrates the importance of tuning spatial scale in liquid-liquid phase separation.

Cortical folding emerges through complex developmental processes, yet what organizes these processes into the characteristic folding patterns of the human brain remains unclear. Here, we developed an imaging-informed whole-brain modeling framework, Grow2Fold, that integrates developmental cortical regionalization, quantitative growth characterization, and mechanics-guided simulation to investigate how cortical growth shapes brain folding patterns. Using 836 fetal and neonatal MRI sessions from 747 participants spanning 21-45 postmenstrual weeks, we identified a reproducible and hierarchical spatiotemporal organization of cortical expansion characterized by coordinated regional growth trajectories. Rapid expansion occurs earlier around peri-central territories and progressively shifts toward frontal, parietal, and temporal association cortices. We encoded these imaging-derived trajectories into region- and age-specific growth laws and incorporated them into anatomically realistic whole-brain mechanical models. The resulting simulations recapitulate major sulcal landmarks and characteristic folding development across global and regional morphological measures. Comparing regional and uniform growth, we find that differential growth drives cortical instability, while heterogeneous developmental growth organizes this instability into anatomically recognizable folding patterns. Combined with subject-specific anatomy, these population-informed growth laws reproduce individualized longitudinal folding trajectories. Together, our results identify spatiotemporally patterned developmental growth as a key spatial constraint on cortical morphogenesis and establish a quantitative bridge from developmental imaging, through mechanical instability, to anatomically organized folding.

Extracellular matrix (ECM)-mediated cell adhesion and tension regulate development and disease through crosstalk with Receptor Tyrosine Kinases (RTKs). RTKs stimulate proliferation, survival, and metabolic reprogramming via Ras/ERK and PI3K/Akt signaling. Emerging evidence identify the biophysical and biochemical integration of integrin-mediated adhesions with RTK signaling, however, the molecular mechanisms regulating this phenotype remain poorly understood. Here, we demonstrate that ECM stiffness potentiates ligand-mediated epidermal growth factor receptor (EGFR) signaling at scales ranging from the single adaptor protein, Grb2 upward to multicellular spheroids and tissues. We show that ECM stiffness-driven RTK signaling is mediated by integrin adhesions that reshape the membrane to induce the formation of RTK signaling hotspots distributed globally across the cell membrane. A parsimonious computational model and biochemical membrane perturbations causally implicate membrane curvature and topography in tension-dependent integrin-EGFR crosstalk. The results demonstrate how the mechanical properties of the ECM direct molecular recruitment and confinement of signaling molecules to regulate growth factor-dependent RTK signaling. Given that many RTK-driven solid tumors develop a stiff stroma, the findings could explain why many therapeutic tyrosine kinase inhibitors develop clinical resistance. The findings suggest clarifying how tension regulates transmembrane receptor signaling could identify unique force-dependent molecular regulators towards which new antitumor therapies could be applied.

Cells use ATP-driven chaperones to extract clients from biomolecular condensates selectively, yet no closed-form criterion links molecular parameters to condensate quality control. We derive a minimal non-equilibrium sorting reactor model that reduces the reaction-diffusion dynamics to one algebraic equation. In the idealized low-load limit, a concentration-dependent partition coefficient produces a saddle-node bifurcation between functional clearance and pathological storage. The critical Damk\"ohler number is \(\mathrm{Da}_{0c}=e^{-2}\) and is independent of the feedback strength. Corrections shift the value but preserve the topology. The model predicts a Kramers-type functional lifetime and recovery costs that scale as \(R^3\). Critical slowing down near the bifurcation provides early-warning signatures for FCS and FRAP.

We investigate the dynamics of the spread of the infection due to a phage- during the range expansion of an E coli colony. We find a rich panoply of spatio-temporal dynamics at large colony-level length scales, arising due to the interplay of mechanical forces, growth and viral infection dynamics at smaller cellular length scales. Strikingly, the phages "surf" the front of the growing bacterial colony by hitchhiking on E coli that are advected due to the growth of bacteria in the colony, resulting in an anisotropic spread of the phages. We identify microscopic processes - of the phage release during lysis (occurring on short millisecond timescales at length scales close to that of a single bacterium) and local nematic alignment of the rod-like E coli bacteria (occurring on timescales comparable to the bacterial growth rate and the length scales of a few bacterial cells) - that enhance the advective effects driving the hitchhiking behavior. Combining our experiments with a mathematical model, we explain the multiple spatio-temporal dynamical regimes - from the coexistence of uninfected bacteria, resistant (lysogenic) cells and phages, to fixed points where the entire population turns resistant or remains uninfected. Our work broaches a new frontier in the investigation of eco-evolutionary processes arising from the mechanical and ecological forces due to the interactions between viruses and their hosts.

Accurate chromosome segregation relies on kinetochores maintaining load-bearing attachments to dynamic microtubule ends, a coupling in which the Dam1 complex ring is central. How the microtubule-end structure and ring-microtubule interactions collectively determine attachment stability and drive force transduction remains unresolved. Here, we use multiscale modeling to show that the ring forms a "fuzzy" complex mediated by a dynamic network of intrinsically disordered regions, enabling both high-affinity binding and free diffusion along the microtubule lattice. Crucially, we find that this prototypic disordered-disordered protein complex provides most of the kinetochore-microtubule stabilization. By contrast, protofilament bending, commonly thought to drive force transduction, is sensitive to the arrangement of protofilaments along the ring and insufficient on its own to establish a robust coupling. The architecture of the dynamic microtubule end unifies these mechanisms by controlling both the progressive loss of fuzzy contacts and the resistance generated by protofilament bending. By accounting for the full conformational ensemble of microtubule structures, our model produces rupture forces similar to those previously measured and explains the distinct, tension-dependent behavior of kinetochore attachments to growing versus shortening microtubule ends. We propose that the Dam1 complex ring acts as a biased-diffusion coupler that probes and gradually remodels the evolving conformational landscape of the microtubule end under tension.